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Abstract. We discuss the maximum kinematical invariance group of the quantum harmonic oscil- 
lator from a viewpoint of the Ermakov-type system. A six parameter family of the square integrablc 
oscillator wave functions, which seems cannot be obtained by the standard separation of variables, 
is presented as an example. The invariance group of the generalized driven harmonic oscillator is 
^ ■ shown to be isomorphic to the corresponding Schrodinger group of the free particle. 

Q" 

Quantum systems with quadratic Hamiltonians are called the generalized harmonic oscillators 
(see, for example, [9], [12], [H], [28], [26], [UJ, [H] and the references therein). These systems 
SpH 1 have attracted substantial attention over the years because of their great importance in many 
^ advanced quantum problems. Examples are coherent states and uncertainty relations, Berry's phase, 
quantization of mechanical systems and Hamiltonian cosmology. More applications include, but are 
not limited to charged particle traps and motion in uniform magnetic fields, molecular spectroscopy 
and polyatomic molecules in varying external fields, crystals through which an electron is passing 
C\J ■ and exciting the oscillator modes, and other mode interactions with external fields. Quadratic 
Hamiltonians have particular applications in quantum electrodynamics because the electromagnetic 
\q \ field can be represented as a set of forced harmonic oscillators |14j . 

in ■ 

The purpose of this paper is to give a simple derivation of the maximum kinematical invariance 
groups of the free particle and harmonic oscillator, which were introduced in Refs. [3], [1], [17], [31] 
and [32] (see also [5] and the references therein), from a unified approach to generalized harmonic 
oscillators (see, for example, [7J, [§], [23] and the references therein). Relations with the correspond- 
ing Riccati and Ermakov-type systems, which seem to be missing in the available literature, are 
emphasized. 



X 



1. Transforming Nonautonomous Schrodinger Equation into Autonomous Form 
Quantum systems described by the one-dimensional time- dependent Schrodinger equation 

where the variable Hamiltonian H = Q (p, x) is an arbitrary quadratic of two operators p = —id/dx 
and x, namely, 

tyt = ~ a (*) ^xx + h (t) x 2 ip - ic (t) xip x -id(t)i)- f (t) xip + ig (t) ip x (1.2) 
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(a, b, c, d, f and g are suitable real- valued functions of time only), are known as the generalized 
(driven) harmonic oscillators. Some examples, a general approach and known elementary solutions 
can be found in Refs. [7], [8], [9], [TO], p2], [H], [26], [37], [46] and [47]. 

In this paper, we shall use the following result established in |23j. 
Lemma 1. The substitution 

e i[a{t)x 2 +8{t)x+K{t)) 

V>= r=77T *^' r )' Z = f3(t)x + e(t), r = - 1 {t) (1.3) 

transforms the non- autonomous and inhomogeneous Schrddinger equation U.S\) into the autonomous 
form 

*Xr = -X« + c^ 2 X (co = 0, 1) (1.4) 

provided that 

dot 

— + b + 2ca + Aaa 2 = c a(3 4 , (1.5) 
dB 

+ (c + Aaa) (3 = 0, (1.6) 



^ 7 +a/3 2 = (1.7) 



and 



Here 



dt 
dS 

— + (c + Aaa) 5 = f + 2ga + 2c a(3 3 e, (1.8) 
dt 

de 

a g = ( g -2a8)P, (1.9) 

^1 = g S -a5 2 + c a/3 2 e 2 . (1.10) 
dt 

I d 

a = . f 1.11 1 

4a n 2a 



(We have changed r — > —t in this paper for convenience. A Mathematica proof of Lemma 1 is 
given by Christoph Koutschan [20]; see also [21]. A special case of the substitution (11. 3p has been 
named the Quantum Arnold Transformation in the recent publications [2], [16] and [27]; see also 
1231. [481 and references therein for the earlier works.) 



The substitution (11.111) reduces the inhomogeneous equation (II. 5p to the second order ordinary 
differential equation 

fj? - t (t) n' + Ao (t) /I = c (2a) 2 /3 4 /i, (1-12) 
that has the familiar time-varying coefficients 

/ n cl' , / \ , ,o d (a 1 d'\ , . 

r{t) = 2c + 4d, a(t) = ab - cd + d 2 + - I -J. (1.13) 

When c = 0, equation (II. 5p is called the Riccati nonlinear differential equation [43], [45] and the 
system ( ll.5p -( fl~TOl) shall be referred to as a Riccati-type system. (Similar terminology is used in 
[39] . [40] for the corresponding parabolic equation.) If Co = 1, equation ( 11.121) can be reduced to a 
generalized version of the Ermakov nonlinear differential equation (see, for example, [9], [13], |25j . 
[41] and the references therein regarding Ermakov's equation) and we shall refer to the corresponding 
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system fll.5p - fll.10p with c 7^ as an Ermakov-type system. Throughout this paper, we use the 
notations from Ref. [23] where a more detailed bibliography on the quadratic systems can be found. 
We have to remind to the reader how to solve the systems (11.5p — fll.lOp (in quadratures) in order to 
make our presentation as self-contained as possible. 

The time- dependent coefficients a , (3 , 7 , 5 , e , k that satisfy the Riccati-type system (|1.5p - 
(JTIUD are given as follows [ZJ, 133, 

"°(*)= (1-14) 
W 4a (t) n (t) 2a (t) ' V ; 

0o (*) = -^|, A (t) = exp (- jf (c (s) - 2d (s)) ds) , (1.15) 
7o(*) = -L^i + ^ d-16) 



2^ (0) /x (t) 2a (0) 



and 



/*0 (*) 



L 



j, < \ d(s) A g (s) 



a (s) J 2a (s) 



ds 

Ms) 



;i.i7) 



2 a(t)A(t) A* a( 3 )<r( 3 )A (a) 



*o(t) = -^^S [t) + S ^^Ms)8 (s))ds (1.18) 



n f t a(s)X(s) (, . d s , A , 
io K)0) V a{s) J 

-0 (t) = 61 (t) - 4 /' a -^§ ( S ) o ( S )) 2 «fa (1.19) 

^0 W ./o (a*o ( s )) 

-2 f 4^- (s) 5 (s)) (f (s) - ^f\g (s)) ds 

Jo M s ) V a (s> J 

(So (0) = — £0 (0) = g (0) / (2a (0)) and k (0) = 0) provided that /x and /i 1 are standard solutions 
of equation f ll,12p with c = corresponding to the initial conditions /i (0) = 0, /i' (0) = 2a (0) 7^ 
and n l (0) 7^ 0, /4 (0) = 0. (Proofs of these facts are outlined in Refs. [7], [H] and [37]. Extensions 
to the nonlinear Schrodinger equations are discussed in Refs. PQ, [6], [38] and [42].) 

Solution of the Riccati-type system in terms of a nonlinear superposition principle is considered 

in 



Lemma 2. The solution of the Riccati-type system / tl.5p -/ TO0j) is given by 

/i(*) = 2/x(0)/x (*)(a(0)+ 7o (*)), (1.20) 

m = - -£§M) = mjWx{t), (1.22) 
^ W 2(a(0)+ 7o (*)) 

/3 2 (0) 

7 = 7 - 77 ^1 (L23) 

4 (a (0) + 7 (t)) 
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and 

(5 t = 5 t — __ 1.24 

2 (a (0) +7 (t)) 

£(t) = (1 .25) 
2 (a (0) +7 (t)) 

«(t) = *{0) + «M- X?% (1-26) 

4 (a (0) +7 (t)) 

m terms o/i/ie fundamental solution \1.1J$ - §1.19) subject to the arbitrary (finite) initial data \i (0) , 
a(0), /3(0) ^0,7(0), 8(0), e(0),k(0). 



The following extension allows to solve the Ermakov-type system 
Lemma 3. TTie solution of the Ermakov-type system U.5\ )- irTTty) w/ien cq = 1 (7^ 0) zs gwen fry 



/i = /i(0)/i ^/3 4 (0) + 4(a(0) + 7o ) 2 , (1-27) 

« = « -/3 2 1 ^ " (0) + 70 5 , (1.28) 

^°/3 4 (0)+4(a(0)+7 ) 2 



= = (,) , (1.29) 



^ (0) + 4 (a (0) + 7o ) 2 Mt) 



and 



1 /3 2 (0) 

7 = 7 (Q)-- a rctan 2 ^ ^ a (0) > (1.30) 



5 _ « e (0) /9 3 (0) + 2 (a (0) + 7o) (j (0) + go) , 31) 

° /5 4 (0) + 4 (a (0) + 7o ) 2 



(1.32) 



2e(0) (a (0)+7o)-/3(0) g (0) + go) 
^/3 4 (0) + 4(q(0)+7o) 2 

« (0) + .0 - e (0) /3 3 (0) , 4 , n ^ ( ° )+ '° , 2 (1-33) 

/3 (0) +4(q(0) + 7 ) 

rn x , ^ 2 (0) /3 2 (0) - (5 (0) + e ) 2 
+ (a (0) + 7 ) 



/3 4 (0) + 4 (q (0) + 7o r 



in terms of the fundamental solution fli. 14\ )- fl.l9\) subject to the arbitrary initial data \i (0) , a (0) , 
0(0) ^0,7(0), 5(0), e(0),«(0). 

Using standard oscillator wave functions for equation ( 11.4ft when Co = 1 results in 

i(a£ 2 +<5:r+K) +i(2ra+l)7 

^ n (x, t) = - == e-^ 2 ' 2 H n (0x + e) , (1.34) 
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where H n (x) are the Hermite polynomials (33] and the solution of the Ermakov-type system (jl.5p - 
fll.lOp is available [23] ■ They are also eigenf unctions, 

(1.35) 



£(*)V„ (x,t) = X(t)[n + -)^ n (x,t), 



of the corresponding dynamic invariant 



E(t) 



X(t) 



(p — 2ax — Sy 

W 2 



+ (fix + ef 



;i.36) 



The Green function of generalized harmonic oscillators 

1 



G(x,y,t) 



exp [i (a (t) x 2 + p (t) xy + 7 (t) y 2 + 5 (t) x + e (t)y + k (t))] 



(1.37) 

has been constructed in Ref. [7]. Evaluation of the Berry phase is discussed in Ref. [36J. 

In the reminder of the paper, we apply these general results to the maximum kinematical invari- 
ance groups of the free Schrodinger equation and of the harmonic oscillator [3T], [32] (see also [5] 
and the references therein). In addition, our approach allows to describe the maximal kinematical 
invariance group of the generalized harmonic oscillators. 



2. Special Cases 



In this section, we show that the maximal kinematical invariance groups of the free Schrodinger 
equation [3T] and the harmonic oscillator [32J (and their isomorphism) can be obtain as special 
cases of the transformation ( 11 .3D . 



2.1. Transformation from the free particle to the free particle. In the simplest case a = 1, 
b = c = d = f = g = co = 0, one finds fi Q = 2t, /i x = 1 and ao = —f3 /2 = 7 = 1/ (4t) , 
= e o — K o — 0. The general solution of the corresponding Riccati-type system is given by 



H (t) = /I (0) (1 + 4a (0)t) 
a(0) 



a(t) 



7(*) = 7(0) 



e(t)=e (0) 



4a (O)t' 

/3 2 (0)t 



l + 4a (O)t' 
2/3 (0)5(0)t 



0{t) 
8(t) 



(3(0) 



1 + 4a (O)t' 

HQ) 

l + 4a (O)t' 



= K (0) 



S 2 (0) t 
1 + 4a (O)t' 



(2.1) 
(2.2) 

(2.3) 

(2.4) 



4a (0) t 

The Ansatz (11. 3p together with these formulas determine the Schrodinger group, namely, the max- 
imum (known) kinematical invariance group of the free Schrodinger equation, as follows [3Tj : 

'a(0)x 2 + 5(0)x-5 2 (0)t 



^(x,t) 



exp 



(0)(l + 4a (0)t) 
p(O)x-20 (0)5(0)t 



1 + 4a (0) t 



+ 5(0) 



1 + 4a (0) t 

/3 2 (0)t 
1 + 4a (0) t 



+ «(0) 



-7(0) 



(2.5) 
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We have established a connection of the Schrodinger group with the Riccati-type system (see also 
[15], [IS], [SO], [22], [11] for the Lie group approach; subgroups of the Schrodinger group and 
their invariants are discussed in [5]). 

The subgroups include the familiar Galilei transformations: 



ip (x, t) = exp 



V V 2 



X {x-Vt + x ,t-t ), (2.6) 



when a (0) = k (0) = 0, /3 (0) = n (0) = 1,7 (0) = t , e (0) = x and 5 (0) = V/2; supplemented by 
dilatations: 

ij(x,t) = X {lx,l 2 t) (2.7) 

with a (0) = 7 (0) = 5 (0) = e (0) = k (0) = 0, \i (0) = 1 and (3 (0) = I] and expansions: 

, , . 1 / mi 2 \ / x £ \ . 
VOM) = r exp z— — -r W -, — (2-8) 



y/l + mt V. 4(l + mt)/ V 1 + mt ' 1 + mt 
(/i(0) = l(^0), //(0)=m), 

- 7^ex P (4)x(-|,-^) (2-9) 

(/i(0) = 0, //(0) = 2#0)) 

with /3 (0) = 1, 5(0) = e (0) = k (0) = . (The symmetry group of the corresponding diffusion 
equation is discussed in [15], [30], [35] and [40J.) 

Here, we derive the Schrodinger group for the free particle as a very special case of Lemma 1. In 
turn, application of (12.51) in (11.31) produces a composition of these two transformations. 

2.2. Transformations from the harmonic oscillator to the free particle. If a = b = 1, 

c = d = f = g = co = 0, the corresponding characteristic equation, p," + 4/x = 0, has two standard 
solutions /i = sin 2t, [i 1 = cos It and 



cos2t a 1 x n 

2 sin 2i sin It 

fi = fi{0) (2a (0) sin 2t + cos 2i) . 

The general solution of the corresponding Riccati-type system takes the form 

2a (0) cos2t - sin2t 



a{t) 
(5 it) 



lit) = 7 (0) 
8(t) - 



2 (2a (0) sin 2t + cos2t)' 

m 

2a (0) sin2t + cos 2t' 

/3 2 (0) sin 2t 



2 (2a (0) sin 2t + cos 2t) 
5(0) 



e(t) = e(0)- 
K (t) = k(0)- 



2a (0) sin2t + cos 2t' 

/3 (0)5(0) sin 2t 



2a (0) sin 2i + cos 2t ' 

5 2 (0)sin2t 
2 (2a (0) sin 2* + cos 2t) 



2.10) 
2.11) 

2.12) 
2.13) 
2.14) 
2.15) 
2.16) 
2.17) 
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subject to arbitrary given initial data. Letting /j (0) = /3 (0) = 1 and a (0) = 7 (0) = S (0) = e (0) 
k (0) = 0, we arrive at the simple substitution 



e -(t/2)x tan2t / tan2t\ 

V x, t = x —57, -5- 2.18 

Vcos2t V cos2 ^ 2 



(see also [E], p], [30], [42J). 



2.3. Transformation from the free particle to the harmonic oscillator. In the simplest case 
a = CQ = l,b = c = d = f = g = 0, the general solution of the corresponding Ermakov-type system 
is given by 



H(t) = n(0) ^4/3 4 (0)t 2 + (4a (0)t+l) 2 , (2.19) 

/3 4 (0)t + a (0) (4a (0) t 
4f3 4 (0) t 2 + (4a (0) t + 1) 



= -q4/ n W2 , /, /nw , ,,2 ■ I 2 ' 20 ) 



0(t) = £ (0) (2.21) 



^4/3 4 (0)t 2 + (4a (0)7+1) 



7 (t) = 7 (0)-iarctan 2 ^ W * ( 2 .22) 
' W V ; 2 4a(0)t + 1 V ; 

= 2e(0)/3 3 (0) t + (0) (4a (0) t + 1) , . 

1 ' 4/3 4 (0)t 2 + (4a(0)t + l) 2 ' 

= e(0) (4a(0)t + l)-2/3 (0) (0) t 

^/J 4 (0)t 2 + (4a (0)t + l) 2 



n(t) = , (0 ) U (4 " (0)t+1) ^ (0)/32(0) -' 2(0) ) (2 25) 

W V ; 4/3 4 (0)t 2 + (4a(0)t + l) 2 V ' ' 

_ e Ae (0) 6 (0) /3 3 (0) 

4£ 4 (0) t 2 + (4a (0)t + l) 2 ' 

When /x (0) = /3 (0) = 1, a (0) = 7 (0) = 5 (0) = e (0) = k (0) = 0, we get the known transformation 
[16], [IT], [IH], [27], [30], 

, . 1 / tx 2 \ ( x 1 \ . 

ip(x,t) = -7- exp j— y . — , - arctan 2t . (2.26) 

^ y ' ' (4t 2 + l) 1/4 V 4t 2 + I J X \y/W+T'2 J 1 ; 

One can easily verify that transformations (12.181) and (12.261) are (local) inverses of each other. 

Equations (11.341) and fl2.19p - fl2.25p provide a six parameter family of the "harmonic oscillator 
states" for the free particle (see also [16] and the references therein). 

2.4. Transformation from the harmonic oscillator to the harmonic oscillator. We consider 
the case a — b — Cq — 1, c = d — f — g — 0. The general solution of the corresponding Ermakov-type 
system is given by 



fi(t) = /i(0) VP (0)sin 2 2t+ (2a(0)sin2t + cos2fT, (2.27) 
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= a (0) cos4t + sin4t (/3 4 (0) + 4a 2 (0) - l) /4 
a ~ /3 4 (0) sin 2 2t + (2a(0)sin2t + cos2f) 2 ' 

0(t) = H ° ] =, (2.29) 



/3 4 (0) sin 2 2t + (2a (0) sin 2t + cos 2ty 

, s 1 /3 2 (0)sin2£ 

7 (t) = 7 (0) - - arctan , N — , 2.30 

' W r V ; 2 2a(0)sin2£ + cos2t' v ; 

5 (0) (2a (0) sin 2t + cos 2t) + e (0) /3 3 (0) sin 2t 

t = J n 9 ) (2.31) 

/3 4 (0) sin 2 2t + {2a (0) sin 2t + cos 2tf 

e (0) (2a (0) sin 2t + cos 2t) - (3 (0) 5 (0) sin 2t 
£ (t) — / 5 (2.32) 

/3 4 (0) sin 2 2t + (2a (0) sin 2t + cos 2i) 2 

« (t ) = M0) + ^2 i £ ' )^ " a ' ^(0)-^(0)^0))-M0)/(0) (233) 

/3 4 (0) sin 2 2t + (2a (0) sin 2t + cos 2t) 2 
+ isin 4i ''(^M- <, (») 



4 /3 4 (0)sin 2 2t + (2a(0)sin2t + cos2t) 2 

(A direct substitution with the help of Mathematica verifies that these expressions indeed satisfy 
the Ermakov-type system.) 

The Ansatz (11. 3p together with the formulas (|2.27|) - (j2.33p explicitly determine the harmonic os- 
cillator group introduced in Ref . [32] . The established connection with the corresponding Ermakov- 
type system allows us to bypass the traditional Lie algebra approach [3], (5], [30], [32]. Combination 
of the transformations (12. 5 p and (12. 9p . (12.181) and ( 12. 26ft implies that the harmonic oscillator group 
is isomorphic to the Schrodinger group of the free particle (see (13. 1 p below) [32]. An analog of the 
similarity transformation (12 .9p takes the form ip (x,t) = %(— x, t — 7r/4) and further identification 
of the corresponding subgroups and the identity element is similar to the case of the free particle. 
We have derived the harmonic oscillator group as a special case of Lemma 1 and one can consider 
a composition of these transformations. 

Use of the explicit solution ( 12 .27[) — (I2.33j) in the wave function (ll.34p results in a six parameter 
(square integrable) family of the "dynamic quantum oscillator states" , which seems cannot be ob- 
tained by the standard separation of variables (the case f3 (0) = 1 and a (0) =7 (0) = 5 (0) = 
s (0) = k (0) = corresponds to the textbook oscillator solution [22], [2S])- The corresponding qua- 
dratic dynamic invariant (I1.36P and the creation and annihilation operators are found (in general) 
in Ref. [36]. A generalization of the coherent states is discussed in Ref. [24] . 

3. The Maximal Kinematical Invariance Group of the Generalized Driven 

Harmonic Oscillators 

The transformation ( 11. 3p from Lemma 1 admits an inversion when the coefficient a (t) does not 
change its sign (see (11.71) for the monotonicity and local time inversion). As a result, the invariance 
group of the generalized driven harmonic oscillator is isomorphic to the Schrodinger group of the 
free particle, 

T = S^TqS, (3.1) 
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thus extending the result of [32] to the corresponding nonautonomous systems (in the classical case, 
see (12.181) and (12.261) for possible operators 5* and S^ 1 and the operator T is defined by (12. 5p and 
( 12. 9p ). The structure of the Schrodinger group of operators To in two-dimensional space-time as a 
semidirect product of SL (2, R) and Weyl W (1) groups is discussed, for example, in Refs. [5], [19] 
and [30]. Further details are left to the reader. 

4. A Conclusion 

Our analysis of the maximal kinematical invariance group of the quantum harmonic oscillator 
provides a six parameter family of solutions, namely fll .34j) and ( I2.2T[) — ( ]2.33j) . for the arbitrary 
initial data of the corresponding Ermakov-type system. These "hidden parameters" disappear af- 
ter evaluation of matrix elements and cannot be observed from the spectrum. How to distinguish 
between these "new dynamic" and the "standard" harmonic oscillator states is thus an open prob- 
lem. (The probability density \ip\ 2 of the solution fll .34j) is moving with time even for the oscillator 
"dynamic ground state" - a simple Mathematica animation reveals such oscillations. This effect, 
quite possibly, can be observed experimentally.) This ambiguity, which is due to the nontrivial 
oscillator maximal kinematical invariance group, then goes further into the coherent states, evalu- 
ation of Berry's phase and dynamic invariants through the established connection with solutions of 
the Ermakov-type system. All of that puts considerations of this paper into a somewhat broader 
mathematical and physical context. 
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